Radar Systems for Glaciology

This chapter deals with radar systems, measurements and instrumentation employed to study the internal core and bedrock of ice sheets in glaciology. The Earth's ice sheets are in Greenland and Antarctica. They cover about 10% of the land surface of the planet. The total accumulated ice comprises 90% of the global fresh water reserve. These ice sheets, associated with the ocean environment, provide a major heat sink which significantly modulates climate. Glaciology studies aim to understand the various process involved in the flow (dynamics), thermodynamics, and long-term behaviour of ice sheets. Studies of large ice masses are conducted in adverse environmental conditions (extreme cold, long periods of darkness). The development of remote sensing techniques have played an important role in obtaining useful results. The most widely used techniques are radar systems, employed since the 1950s in response to a need to provide a rapid and accurate method of measuring ice thickness. Year by year, polar research has become increasingly important because of global warming. Moreover, the discovery of numerous subglacial lake areas (water entrapped beneath the ice sheets) has attracted scientific interest in the possible existence of water circulation between lakes or beneath the ice (Kapitsa et al., 2006; Wingham et al., 2006; Bell et al., 2007). Recent studies in radar signal shape and amplitude could provide evidence of water circulation below the ice (Carter 2007, Oswald and Gogineni 2008). In this chapter the radar systems employed in glaciology, radio echo sounding (RES), are briefly described with some interesting results. RES are active remote sensing systems that utilize electromagnetic waves that penetrate the ice. They are used to obtain information about the electromagnetic properties of different interfaces (for example rock-ice, ice-water, seawater-ice) that reflect the incoming signal back to the radar. RES systems are characterized by a high energy (peak power from 10 W to 10 KW) variable transmitted pulse width (about from 0.5 ns to several microseconds) in order to investigate bedrock characteristics even in the thickest zones of the ice sheets (4755 m is the deepest ice thickness measured in Antarctica using a RES system). Changing the pulse length or the transmitted signal frequencies it is possible to investigate particular ice sheet details with different resolution. Long pulses allows transmission of higher power than short pulses, penetrating the thickest parts of the ice sheets but, as a consequence, resolution decreases. For example, the GPR system, commonly used in geophysics for rock, soil, ice, fresh water, pavement and structure characterization, employs a very short transmitted pulse (0.5 ns to 10 ns) that allow detailing of the shallow parts of an ice sheet (100-200 m in depth) (Reynolds 1997). Consequently, in recent years, GPR systems are also employed by explorers to find hidden crevasses on glaciers for safety. RES surveys have been widely employed in Antarctic ice sheet exploration and they are still an indispensable tool for mapping bedrock morphologies and properties of the last unexplored continent on Earth. The advantage of using these remote sensing techniques is that they allow large areas to be covered, in good detail and in short times using platforms like aeroplanes and surface vehicles

Geothermal flux and basal melt rate in the Dome C region inferred from radar reflectivity and heat modelling

Abstract. Basal melt rate is the most important physical quantity to be evaluated when looking for an old-ice drilling site, and it depends to a great extent on the geothermal flux (GF), which is poorly known under the East Antarctic ice sheet. Given that wet bedrock has higher reflectivity than dry bedrock, the wetness of the ice–bed interface can be assessed using radar echoes from the bedrock. But, since basal conditions depend on heat transfer forced by climate but lagged by the thick ice, the basal ice may currently be frozen whereas in the past it was generally melting. For that reason, the risk of bias between present and past conditions has to be evaluated. The objective of this study is to assess which locations in the Dome C area could have been protected from basal melting at any time in the past, which requires evaluating GF. We used an inverse approach to retrieve GF from radar-inferred distribution of wet and dry beds. A 1-D heat model is run over the last 800 ka to constrain the value of GF by assessing a critical ice thickness, i.e. the minimum ice thickness that would allow the present local distribution of basal melting. A regional map of the GF was then inferred over a 80 km  ×  130 km area, with a N–S gradient and with values ranging from 48 to 60 mW m−2. The forward model was then emulated by a polynomial function to compute a time-averaged value of the spatially variable basal melt rate over the region. Three main subregions appear to be free of basal melting, two because of a thin overlying ice and one, north of Dome C, because of a low GF

Spatial and temporal distributions of surface mass balance between Concordia and Vostok stations, Antarctica, from combined radar and ice core data: first results and detailed error analysis

Results from ground-penetrating radar (GPR) measurements and shallow ice cores carried out during a scientific traverse between Dome Concordia (DC) and Vostok stations are presented in order to infer both spatial and temporal characteristics of snow accumulation over the East Antarctic Plateau. Spatially continuous accumulation rates along the traverse are computed from the identification of three equally spaced radar reflections spanning about the last 600 years. Accurate dating of these internal reflection horizons (IRHs) is obtained from a depth-age relationship derived from volcanic horizons and bomb testing fallouts on a DC ice core and shows a very good consistency when tested against extra ice cores drilled along the radar profile. Accumulation rates are then inferred by accounting for density profiles down to each IRH. For the latter purpose, a careful error analysis showed that using a single and more accurate density profile along a DC core provided more reliable results than trying to include the potential spatial variability in density from extra (but less accurate) ice cores distributed along the profile. The most striking feature is an accumulation pattern that remains constant through time with persistent gradients such as a marked decrease from 26 mm w.e. yr(-1) at DC to 20 mm w.e. yr(-1) at the south-west end of the profile over the last 234 years on average (with a similar decrease from 25 to 19 mm w.e. yr(-1) over the last 592 years). As for the time dependency, despite an overall consistency with similar measurements carried out along the main East Antarctic divides, interpreting possible trends remains difficult. Indeed, error bars in our measurements are still too large to unambiguously infer an apparent time increase in accumulation rate. For the proposed absolute values, maximum margins of error are in the range 4 mm w.e. yr(-1) (last 234 years) to 2 mm w.e. yr(-1) (last 592 years), a decrease with depth mainly resulting from the time-averaging when computing accumulation rates

Fast geophysical prospecting applied to archaeology: results at «Villa ai Cavallacci» (Albano Laziale, Rome) site

The present essay is the result of a cooperative work between geophysicists and archaeologists in which the authors carried out an integrated geophysical prospecting in an archaeological site near Rome. This paper describes the methodology and the results of a geophysical survey carried out on Villa ai Cavallacci, an ancient roman building in Albano Laziale (Rome) discovered in the late seventies. It is often possible to obtain very important results planning a fast geophysical survey opportunely; within this framework (due to the fact that an archaeological excavation was planned in a short time), an integrated geophysical techniques survey (GPR, magnetic, and geoelectric tomography) has been carried out on the areas indicated by the archaeologists. Even if the described geophysical survey should be considered only a first step analysis, the data pointed out some very interesting features confirmed by the excavation

Location of a new ice core site at Talos Dome (East Antarctica)

In the frame of glaciology and palaeoclimate research, Talos Dome (72°48lS; 159°06lE), an ice dome on the East Antarctic plateau, represents the new selected site for a new deep ice core drilling. The increasing interest in this re- gion is due to the fact that the ice accumulation is higher here than in other domes in East Antarctica. A new deep drilling in this site could give important information about the climate changes near the coast. Previous papers showed that the dome summit is situated above a sloped bedrock. A new position on a relatively flat bedrock 5-6 km far from here in the SE direction was defined as a possible new ice core site for an European (Italy, France, Swiss and United Kingdom) drilling project named as TALDICE (TALos Dome Ice Core Project). This point, named as ID1 (159°11l00mE; 72°49l40mS), became the centre of the Radio Echo Sounding (RES) flight plan during the 2003 Italian Antarctic expedition, with the aim of confirming the new drilling site choice. In this paper 2001 and 2003 RES data sets have been used to draw a better resolution of ice thickness, bottom morphology and internal layering of a restricted area around the dome. Based on the final results, point ID1 has been confirmed as the new coring site. Fi- nally, the preliminary operations about the installation of the summer ice core camp (TALDICE) at ID1 site carried out during the XX Italian Antarctic expedition (November 2004-December 2005) are briefly described

Active and capable fault? The case study of Prata D'Ansidonia (L'Aquila, Central Apennine)

The study deals with the morphogenetic meaning of several linear scarps that carved the paleo-landsurface of Valle Daria, an extended geomorphological feature located between Barisciano (AQ) and Prata D'Ansidonia (AQ). These villages are situated in the southern termination of the L'Aquila intermontane basin (one of the largest basin of the central Apennines), nearby the epicentral area of the 6th April 2009 earthquake (Mw 6.1). These scarps, up to 3 meters high and up to 1.5 km long, define narrow/elongated flat-bottom depressions, filled by colluvial deposits. These depressions are carved into fluvial-deltaical conglomerates, dated back to the lower Pleistocene. Even if different authors have interpreted these shapes as a paleodrainage or secondary faults, a morphometrical study of the Valle Daria paleo-landsurface provided several information which cast doubt on these two interpretations. In order to better understand the nature and the state of activity of these lineaments, geological, geomorphological and geophysical surveys were carried out. A paleoseismological trench pointed out two events of deformation. The curvilinear shape of the shear plane seems to be related to a slow deformation, attributable to collapse-phenomena. Three GPR profiles, two ERT profiles and two microgravimetrical profiles seem to corroborate this interpretation. Therefore, this study permits to attribute the genesis of these scarps to tectono-karstic phenomena, excluding the presence of an active and capable fault.Published346-3494T. Sismologia, geofisica e geologia per l'ingegneria sismicaN/A or not JC

Comparison of measurements from different radio-echo sounding systems and synchronization with the ice core at Dome C, Antarctica

We present a compilation of radio-echo sounding (RES) measurements of five radar systems (AWI, BAS, CReSIS, INGV and UTIG) around the EPICA Dome C (EDC) drill site, East Antarctica. The aim of our study is to investigate the differences of the various systems in their resolution of internal reflection horizons (IRHs) and bedrock topography, penetration depth, and quality of imaging the basal layer. We address the questions of the compatibility of existing radar data for common interpretation, and the suitability of the individual systems for Oldest Ice reconnaissance surveys. We find that the most distinct IRHs and IRH patterns can be identified and transferred between most data sets. Considerable differences between the RES systems exist in range resolution and depiction of the basal layer. Considering both aspects, which we judge as crucial factors in the search for old ice, the CReSIS and the UTIG systems are the most valuable ones. In addition to the RES data set comparison we calculate a synthetic radar trace from EDC density and conductivity profiles. We identify ten common IRHs in the measured RES data and the synthetic trace. The reflection-causing conductivity sections are determined by sensitivity studies with the synthetic trace. In this way, we accomplish an accurate two-way travel time to depth conversion for the reflectors, without having to use a precise velocity-depth function that would accumulate depth uncertainties with increasing depth. The identified IRHs are assigned with the AICC2012 time scale age. Due to the isochronous character of these conductivity-caused IRHs, they are a means to extend the Dome C age structure by tracing the IRHs along the RES profiles

Radio-echo sounding at Dome C, East Antarctica: A comparison of measured and modeled data

The internal layering architecture of ice sheets, detected with radio-echo sounding (RES), contains clues to past ice-flow dynamics and mass balance. A common way of relating the recorded travel time of RES reflections to depth is by integrating a wave-speed distribution. This results in an increasing absolute error with depth. We present a synchronization of RES-internal layers of different radar systems (Alfred Wegener Institute, Center for Remote Sensing of Ice Sheets, Istituto Nazionale di Geofisica e Vulcanologia, British Antarctic Survey and University of Texas Institute for Geophysics) with ice-core records from the Antarctic deep drill site Dome C. Synthetic radar traces are obtained from measurements of ice-core density and conductivity with a 1D model of Maxwell’s equations. The reflection peaks of the different radar systems’ measurements are shifted by a wigglematching algorithm, so they match the synthetic trace. In this way, we matched pronounced internal reflections in the RES data to conductivity peaks with considerably smaller depth uncertainties, and assigned them with the ice-core age. We examine the differences in shifts and resolution of the different RES data to address the question of their comparability and combined analysis for an extensive age-depth distribution